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Abstract 


In this study, we propose a family of iterative procedures with no deriva- 
tives for calculating multiple roots of one-variable nonlinear equations. We 
also present an iterative technique to approximate the multiplicity of the 
roots. The new class is optimal since it fits the Kung—Traub hypothesis and 
has second-order convergence. Derivative-free methods for calculating mul- 


tiple roots are rarely found in literature, especially in the case of one-step 
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methods, which are the simplest ones in terms of their structure. Moreover, 
this new family contains almost all the existing single-step derivative-free 
iterative schemes as its special cases, with an additional degree of freedom. 
Several results are used to confirm its theoretical order of convergence. 
Through the complex discrete dynamics analysis, the stability of the sug- 
gested class is illustrated, and the most stable methods are found. Several 
test problems are included to check the performance of the proposed meth- 
ods, whether the multiplicity of the roots is estimated or known, comparing 


the numerical results with those obtained by other methods. 
AMS subject classifications (2020): Primary 65H05; Secondary 37D99. 


Keywords: Nonlinear equations; Derivative—free iterative method; Multiple 


roots; Order of convergence; Stability. 


1 Introduction 


One well-known derivative-free scheme for obtaining a multiple root x* of 
a nonlinear complex-valued function f(x), « € C is the modified Traub-— 


Steffensen’s approach [19]. The structure of this scheme is as follows: 


f(t) 


x =%—-—m ,» €20,1,23..., 1 
“a ‘ f ln, 24] ” 


| = Lada flan) 


Nte— Xt 


the first-order divided difference at point x; and m = a; + Bf(a1), 8 4 90, 


B €R. For m > 2, the order of convergence of this scheme is quadratic. 


where m denotes the known multiplicity of «* and f[n, xv: is 


Scheme (1) is the derivative-free version of the familiar modified Newton’s 
method [15] as it is obtained from this method by estimating the derivative 
f' (xz) by fm, x4]. Centered on these two schemes, many methods [6, 13, 
16, 18, 22, 11, 17, 21, 7, 10] have been derived in earlier studies. Among 
these, the second-order schemes presented by Kumar, Sharma, and Argyros 
[9] and Kansal et al. [8], are the nicest ones as they possess the following 
features: (i) the schemes are derivative-free, (ii) they are optimal in terms of 
the Kung-Traub conjecture (see [10]), and (iii) the schemes are single-step, 


which make them easy to implement. The general type of optimal derivative- 
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free procedure for multiple roots defined by Kumar, Sharma, and Argyros [9] 
is given by 
Tt4+1 = Le — H(®;), (2) 


f(t) 
fle, ve] 


main aspect of this scheme is that it holds Traub-Steffensen’s method (1) as 


where ®; = and H(®,;) is a single variable weight function. Other 


a special case. 


The scheme from Kansal et al. [8] has the following structure: 


— fb OF (mn) + af (a1) 
Tt4+1 = Ut fim. 2] ’ (3) 


where a € R is a free parameter. For a = 1, this scheme reduces to Traub-— 
Stefensen’s method (1). 

Motivated by these multiple root-finding schemes, we aim to design a new 
derivative-free iterative family that is more general than class (2). Therefore, 
the new class is a single-step optimal scheme with second-order convergence, 
and it has an additional degree of freedom. In addition, we include an itera- 
tive process that allows us to estimate the multiplicity when it is not known. 

In order to choose the most stable members of the proposed class, we now 
introduce some dynamic terms used in this manuscript (see, for instance, 
[2]). To achieve this aim, the key concept is the rational function. Usually, a 
rational function T' is found by applying the proposed class on a low-degree 
polynomial p(z) with multiple roots. The properties of this rational function, 
the existence of fixed points or periodic orbits, and their repulse or attracting 
character, among others, allow us to detect the best elements of the class or 
to compare known schemes (see, for example, [4, 5]). 

Let us consider a rational function T : C —> GC, being C the Riemann 


sphere. Then, we define the orbit of a point z as follows: 
{z,T(2),T?(2),-.., T(z), 


where T* denotes the kth composition of T with itself. 

We analyze T by characterizing the starting points from the asymptotical 
performance of their orbits. In these terms, a point z is a fixed point of T 
if T(z) = z; it is a p-periodic point, p > 1, if T?(z) = z and T*(z) ¥ z for 
k<p. 
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The basin of an attractor € is defined as 
A(é) = {z € ©: T"(z) > En > oo}. (4) 


The Fatou set of T, is the set of points whose orbits tend to an attractor 
(fixed point, periodic point, or infinity). It’s complementary in C is the Julia 
set. So, the basin of attraction of any fixed point belongs to the Fatou set 


and the boundaries of these basins of attraction belong to the Julia set. 


On the other hand, a point z is called a critical point of T if T’(z) = 0. 
The asymptotic behavior of the critical points is a key fact for analyzing 
the stability of the method: A classical result from Fatou and Julia (see, for 
instance, [1]). Each immediate basin of attraction holds at least one critical 
point; that is, in the connected component of the basin of attraction holding 


the attractor, there is also a critical point. 


Moreover, a fixed point of T, z, is called attracting if |T’(z)| < 1, or 
super attracting if |T’(z)| = 0; it is repulsive if |T’(z)| > 1; and parabolic if 
|T’(z)| =1. 

Let us also remark that when fixed and critical points are not equivalent to 
the roots of p(z), then they are called strange fixed and free critical points, 
respectively. The asymptotical performance of free critical points is a key 
element in this study, as when they do not belong to the basin of attraction 
of the searched zeros. Then, they lay in the basin of an attracting element 


(a periodic orbit or a strange fixed point) that should be avoided in practice. 


Indeed, when T also depends on one or several parameters r;, 7 = 
1,2,...,k, and the stability of a fixed point z is analyzed, |T’(z,11,...,7Tx)| 
is not a scalar but a function of r; for 7 = 1,2,...,k. In this case, 
|T"(z,T1,.--,7k)| is called a stability function of the fixed point, and it gives 


us the character of the fixed point in terms of the value of r;, i =1,2,...,k. 


Let us summarize the rest of the paper: Section 2 includes the construc- 
tion and convergence analysis of the new family. Also, in Section 2, special 
cases of the proposed class are found. To examine the stability of the new 
methods, we present the analysis of the discrete dynamical system associated 
with the different sub-classes defined in Section 3, selecting the most stable 


ones. Section 4 is devoted to the application of proposed methods on a num- 
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ber of numerical problems to back up the theoretical conclusions. In Section 


5, we conclude with some final observations. 


2 Establishment of second-order class 


Let us define a new derivative-free Traub—Steffensen-type class of iterative 


schemes for finding multiple roots as 


Le SS fe L(®:), (5) 


where L(Q;) is a single-variable weight function, 0; = mos VER, 
and m = 21+ Pf(t), B40, BER. 

We show in the following result that the proposed scheme (5) achieves 
maximum second-order convergence for every 8 4 0, without using any 
derivative evaluation. Moreover, we show that the conditions to be imposed 


on the weight function are directly related to the multiplicity of the root. 


Theorem 1. Let f : DC C-—C be an analytic function in the open set 
D, that holds a multiple root «*, with multiplicity m. Then, the class of 


iterative schemes (5) has second-order of convergence, if 
L(0)=0, L'(0)=m, |L"(0)| < ov, (6) 
and therefore it satisfies the error equation 


1 
etl = 53 (2m(y + Bi) — L""(0)) ef + O(e?), 


1 FO) 
m+1 fl )(a2r) ° 


where e; = 2; — x* and By = 


Proof. Along Taylor’s series expansion of f(a) and f(7;) around «* that 
has multiplicity m, it is assumed that f(a*) = f’(a*) = f"(a*) =--- = 
fo"-Y (a*) = 0 and f(™(2*) £0. Also, let us denote B; by 
Be m! fore (a*) 
“ (m+i)l f(a) ’ 


for 1 = 1,2,3,.... Then, 
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(m) (4* 
f (at) = Pete) on (1 + Byey + Bye? + B3e? + ote!)). (7) 
(m) (4% 
f(m) = Pete) om (1 + Bren, Boe;, B3e}, ole), (8) 
where 
Ff (2*) 


en = Nt — g=et+ B Ee,” (1 + Byex + Boe? + B3e? + ote?)) , 


m! 
By replacing (7)—(8) in the definition of ©,, we obtain 


i 
Bee Sy Be Ole), (9) 


m  m? 


From (9), it can be inferred that ©; tends to zero when ¢ tends to infinity. 


Then, by using MacLaurin’s theorem on the weight function [(©;), we obtain 
1 
L(®,) = L(0) + L’(0)©, + al (0) 8% + 0(03). (10) 
By imposing (9) and (10) in (5), we get 


L (0) 


m 


En Fie (: )e ! = (2(7+ Bi) LE (0)—L°) (0) )e? +0(€?). 
(11) 
From this error equation, we deduce that schemes belonging to class (5) 
possess, at least, second-order of convergence, provided that L(0) = 0 and 
1— £© 0, that is, L’/(0) =m. 
Finally, we find the error equation 


1 
ern = 5g (2m(7 + Br) - L(0)) e? + OC), (12) 


provided |L/’(0)| < oo. Hence, the proposed class (5) converges to the m- 


multiple root x* with second-order convergence. 


Remark 1. It is known that E = pe is the computational efficiency in 
which p is the order of convergence of the method under consideration and w 
is the number of functional evaluations per iteration; see [14]. The E-value 
corresponding to the newly proposed scheme is = 22 = 1.414, which is the 
same as that of modified Newton’s method and modified Traub—Steffensen’s 
method. 
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Remark 2. For 7 = 0, the proposed class (5) reduces to scheme (2), and 
in addition, if we also have L(O,) = mOz, then we get the modified Traub- 
Steffensen’s method (1). 


2.1 Special cases 


In this section, we show that the number of special cases for the new presented 
family is equal to the particular forms of the function L(@;,), fulfilling the 


convergence conditions of Theorem 1. A few of them are listed in Table 1. 


Table 1: Some special cases of class (5) 


Method Weight function Subcases 

TM-1 = L(©) = mO,(1+a0,) a= 4h, tey1 = 1% — MO,(1 + YE) 
TM-2 = L() = ~@& a= HF, G41 = 0 — 328e 

TM-3 L(®,) = _@Ss a= >, ti = — POS 
TM4 L(@)=qp@@gn~ = ote = — GA 
TM-5 (0) = meit90 a =k, wey = a — Cd 


3 Dynamical analysis 


In order to arrange this analysis, we apply our proposed sub-families TM-1 
to TM-5 on the nonlinear function p(z) = (z — 1)™(z+ 1), with z=lasa 
multiple root of multiplicity m and a simple root at z = —1. In all cases, 
y=land f= 10: qualitatively, similar results are found for other values 
of y and £, but the values of 8 close to zero make the schemes more stable, 
due to the estimation of the derivative provided by the first-order divided 
difference is better. This is the most simple nonlinear function containing 
two roots, one simple and one m-multiple. Although the results cannot be 
directly extrapolated to any nonlinear function, several analyses on different 


nonlinear problems confirm, in the numerical section, these results. 
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In what follows, we analyze the rational functions related to the particular 
cases of (5) presented: TM-1 to TM-5. 


3.1 Fixed points and stability of rational TM cases 


When the classes TM-1, TM-2, TM-3, TM-4, and TM-5 are applied on the 

polynomial p(z) with multiplicity m = 2, rational functions T;(z,a), i = 

1,2,...,5 are, respectively, obtained. They all depend on the parameter of 

the class a, considered a free complex value. These rational operators are 
Ni(z,a) 

~ (25 — 24 4 2823 + 922? 4 271z +9) 


being Ni(z,a) = —20000a (z? — 1)” + 241 — 2210 4.5729 + 12828 + 94227 + 
48282 + 1822225 + 3176824 + 2649728 + 2147822 + 542812 + 1800, 


Tol, 0) — 4254+ 28)2" + (9 100a)z + 28 ~ 25 + 2824 + 71a” + 200 
Zz = 
pee 100a (z?2 — 1) + 25 — 24 4 2823 + 922? + 271249 , 
N: 
T3(z, a) = 3(,4) 


(25 — 24 42823 4+ 9222 + 271249)’ 


where N3(z,a) = —40000a (2? — 1)? + 24! — 2210 + 5729 + 12828 + 94227 + 
48282° + 182222° + 3176824 + 2649729 + 214782? + 542812 + 1800, 


Na(z, a) 


T4(z, a) = 2? 
(100V/2a (x? — 1) +25 — a4 + 28x3 + 92x? + 271x + 9) 


being N4(z,a) = 21! — 2219 + 5729 + 8(25/2a + 16) a8 + (942 — 200/2a) ax" + 
(5400V2a+4828)a° +2(10000a?+9300V2a+9111)a°+8(6075V2a+3971)a4+ 
(—40000a? — 16600./2a + 26497) + (21478 — 54200\/2a)a? + (20000a? — 
1800\/2a + 54281) + 1800, and 


N3(z,@ 


Ts(z,a) = Ds Ds(z, a) a 


? 


) 

) 
where N5(z,a) = —18200+ (54281 — 1800a)z+ (61478 — 54200a)z? + (26497 — 
16600a) 23 +8(1471+6075a)z4 +.6(3037+3100a)z° + (4828 +5400a)z6 + (942 — 
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200a)z? +8(16 + 25a)28 +5729 — 2219 +.2" and Ds(z,a) = (9+271z+9227 + 
2823 — 24 + 25)(9 + 2712 + 9227 + 2823 — z4 + 2° + 200a(—1 + z7)). 

By solving each equation T;(z,a) = z, i = 1,2,...,5, the fixed points 


of the ith rational function are obtained. Two of them are the roots of 
polynomial p(z); the rest of them, if they exist, are strange fixed points. 
The asymptotical behavior of all the fixed points (both multiple and simple, 
strange or not) plays a key role in the stability of the iterative methods 
involved, as the convergence to fixed points different from the roots means 
an important drawback for an iterative method; so, we proceed below with 
this analysis. 

In order to study the stability of the fixed points z" of T;(z,a), we cal- 
culate its stability function, It ce a)|. This function gives us information 
about the asymptotical behavior of the point in terms of the value of a. In 
general, the stability of other fixed points than the multiple roots of p(z) 
depends on the value of parameter a. From these rational functions, the 


following result can be stated. 


Theorem 2. Rational function T;(z,a) has z = 1 and z = —1 as super- 
attracting and parabolic fixed points, respectively. There can also exist 
strange fixed points, denoted by ri (a), where the index i = 1,2,...,5 corre- 
sponds to the rational operator T;(z,a) and the index j = 1,2,...,5 corre- 
sponds to the fixed point as the root of fifth-degree polynomial q;(t) described 


below: 


i = 1: Roots of polynomial qi (t) = 9—100a+271t+(92+100a)t? +28¢3—t*+t° 
are the strange fixed points of T\(z, a). 


i = 2: There are no strange fixed points of T(z, a). 


i = 3: When 73(z,a) is analyzed, the strange fixed points are the roots of 
q3(t) = 9 — 200a + 271¢ + (92 + 200a)t? + 28t3 — 4 + 2°. 


i= 4: T,(z,a) has as strange fixed points the roots of q4(t) = 9 + 271¢ + 
92¢? + 2843 — 44 2°. 


i =5: Strange fixed points of Ts(z,a) are the roots of g5(t) = —91 + 271t+ 
1922? + 28t? — ¢4+1°. 
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Proof. By solving the equation T;(z,a) = z, i = 1,2,3,5, it is found that, in 


the case when 7 = 1, 


Ti(z,a) =z Ny(z,a) = 2 (2? — 24 + 2827 + 922? + 2712 4 9)”, 
or, equivalently, 
(2? —1)(9+2712+92z7 +2823 —z*+z°+100a1(z?—1)) = 0 © (z?—1)q1(z) = 
For 7 = 2, it can be checked that 


T(z,a) =z26 24 -1=0. 


In a similar way, the statements for 7 = 3, 7 = 4, and 7 = 5 are obtained. 


So, it has been proven that the multiple root is always a superattracting 
fixed point, whereas the simple root is parabolic. Hence, it always lays in 
the Julia set and would never be reached exactly by the iterative method 
(although it can be very close). For i = 1,3,4,5, there exist strange fixed 
points of their respective operators, T;(z,a). Regarding the stability of these 


fixed points, the following remark can be checked. 


Remark 3. Operators T;(z,a) for i = 1,3,4,5, have attracting strange fixed 


point at any complex value of a, except in a disk |a| < r with 
e r= 1.12 for T,(z, a); 
e r 0.56 for T3(z, a); 


e r 0.75 for Ty(z, a); 


e r 0.11 for Ts(z, a). 


These small regions, where all the strange fixed points are repulsive, can 
be observed in Figure 1. In it, the code color is as follows: One point of the 
mesh is represented in orange color if there exists at least one strange fixed 
point that is attracting this value of the parameter. So, if one point of the 
mesh is not colored, it means that for this value of a, all the strange fixed 


points are repulsive. 
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(a) IT{(r7, 4)| (b) |T3(r3,a)| (c) ITi(r4, a) | (d) |T4(r3,a)| 


Figure 1: Union of stability functions for rational operators on strange fixed points rl (a), 
4=1,3,4,5 and j =1,2,...,5. 


3.2 Dynamical planes 


Each value of the parameter corresponds with a member of each subclass 
of iterative methods. When a is fixed in a rational function T;(z,a), 
i = 1,2,...,5, the performance of the iterative process can be visualized 
in a dynamical plane. It is obtained by iterating the chosen element of the 
family under study and by using each point of the complex plane as an initial 
estimation. In this section, we have used a mesh of 1000 x 1000 points. 

We represent in blue color those points whose orbit converges to infinity, 
in orange the points converging to the multiple root z = 1 (with a tolerance 
of 10~%), and in other colors (green, red, etc.) those points whose orbit 
converges to one of the fixed points (all fixed points appear marked as a 
white star in the figures if they are attractive or by a white circle if they are 
repulsive or parabolic). Moreover, a point appears in black if it reaches the 
maximum number of 500 iterations without converging to any of the fixed 
points. The routines used appear in [3]. 

Firstly, we plot the dynamical planes related to rational functions T;(z, a), 
i = 1,2,3,5 for a = 0.05. These planes can be visualized in Figure 2. From 
Theorem 2 and Remark 3, only the multiple root of p(z) is an attracting fixed 
point. 

From the dynamical planes of Figure 2, it is again deduced that the only 
basin of attraction is that of the multiple roots; moreover, if there exist 
strange fixed points, they lay in the Julia set, as per their repulsive character. 


Black areas correspond to divergent performance. The wideness of the basins 
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(a) T(z, 0.05) (b) T(z, 0.05) (c) T3(z, 0.05) 


(d) T(z, 0.05) (e) T5(z, 0.05) 


Figure 2: Dynamical planes corresponding to rational operators under study for a = 0.05. 


of attraction composing the Fatou set is similar in all cases but slightly better 
in the case of T5(z,0.05). 

In Figure 3, the dynamical planes corresponding to a = 5 are presented. 
For this value of the parameter, only T> operator does not attract strange 
fixed points, from Theorem 2 and Remark 3. 

We observe in Figure 3e that T5(z,5) has three attracting strange fixed 
points, with their corresponding basins of attraction in green, red, and dark 
aubergine color, respectively. Also, T;(z,5), Ts(z,5), and T3(z,5) show green 
basins of attraction of strange fixed points. The latter case is especially 
inefficient, as the basin of attraction of the multiple roots is really small. 
The basin of attraction of the multiple roots in the case of the operator 
T2(z,5) is the only one in the dynamical plane, although it is narrower than 
in the case of a = 0.05. 

Figure 4 corresponds to the dynamical planes of rational functions T;(z, —1), 


i= 1,2,...,5. Cases i = 3, 1 = 4, and i = 5 show basins of attraction of 


Iran. J. Numer. Anal. Optim., Vol. 14, No. 1, 2024, pp 291-314 


303 Title Suppressed Due to Excessive Length 


(a) Ti (z, 5) (b) T(z, 5) (c) T3(z, 5) 


(d) T4(z, 5) (e) Ts (z, 5) 


Figure 3: Dynamical planes corresponding to rational operators under study for a = 5. 


strange fixed points. Meanwhile cases i = 1 and 7 = 2 only converge to the 
multiple roots, being the last one the widest. 

We conclude from the analysis made that the most subclass is T(z, a), 
although the rest of the operators perform adequately for small absolute 


values of a. 


4 Numerical results 


In this section, we use numerical examples to confirm the performance of 
the newly created iterative approaches as well as to validate the theoretical 
results presented in the previous sections. Furthermore, we compare our 
new proposed methods TM-1-TM-5 with the second-order methods given 
by equations (1) and (3). We denote these methods by MNM and KM, 
respectively. For the scheme (3), we choose the parameter a = 6/7. All the 
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(a) Ti (z, —1) (b) Ta(z, =1) (c) T3(z, —1) 


(d) Ty(z, -1) (e) T5(z,—-1) 


Figure 4: Dynamical planes corresponding to rational operators under study for a = —1. 


computations for the above-mentioned methods are performed using software 
Mathematica with 4096 significant digits as multiple precision arithmetic. In 


numerical results, we have calculated the following terms: 


e The number of iterations ¢ required to obtain the desired root fulfilling 


the condition that |a441 — 2;| + |f(a41)| < 107-°°. 
¢ Estimated error e441 = |2441 — 2;| in the last three iterations. 
e Residual error for the considered function, that is, | f(v1+41)|. 


e The computational order of convergence COC [20], which is given by 


t=1,2,.... (13) 


The problems used for validating the theoretical results and comparison are 


mentioned in Table 2. The numerical results obtained for these problems 
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have been presented in Tables 3-6. From them, we see that the new proposed 


methods are, at least, as stable as those existing methods used for comparison. 


Also, the presented results have verified the theoretical results proven so far. 
Moreover, from Tables 3-6, it is observed that methods TM-2, TM-4, and 


TM-5 have performed rightly well in all the problems. 


Furthermore, the 


performance of the methods TM-1 and TM-3 is either at par or better than 


the results obtained for the existing methods. 


Table 2: Test functions considered, their corresponding roots and multiplicity 


Test Functions Xo Root m 
fi(x) = x3 — 5.222? + 9.08252 — 5.2675 1.8 1.75000... 2 
2 
2 

fo(x) = (tan-X(2 ~l) +e” — 3) 1 1.05653... 2 

f3(x) = (sinxcosx — x3 + 1)9 14 1.11707... 9 

fa(x) = ((x — 1)? — 1)*° 2.1 2.00000... 50 

Table 3: Numerical results for f1 (a) 

Methods t ler—1| lez| ler4a| lf (te41)| COC 
TM-1(y=0) 8 5.45x 10715 4.93x 10778 4.04x 10754 2.21 x 107-743 2.00000 
(y=-1 8 145x107 3.38 x 10-79 1.84x« 107° 8.94 x 10-728 2.00000 
TM-2 (y=0) 8 3.25x107'5 1.73x 10778 4.92x107°° 4.76 x 107-217 2.00000 
(y=-1 8 825x107 1.08 x 10-29 1.87x 10757 9.27 x 10-727 2.00000 
TM-3 (y=0) 8 1.66x 10715 4.46x 10779 3.22x 107° 8.40 x 10-27? 2.00000 
(y=-1 8 3.97x 10-16 2.47x 10-89 9.55 x 10-59 6.11 x 10-78? 2.00000 
TM-4 (y=0) 8 9.41x1071® 1.41 x 10779 3.19x 10757 7.90 x 10-276 2.00000 
(y=-1 8 2.13x 107-1 7.02 x 10-* 7.63 x 107-8? = 2.42 x 10778 2.00000 
TM-5 (y=0) 8 9.76x 10716 1.52x 10-79 3.69x 107-57 1.42 x 107-2? 2.00000 
(y=-1 8 222x107 7.62 10-%! 8.98 x 107° 4.67 x 10-78 — 2.00000 
MNM 8 618x107 637x 10-7 6.76x 10-4 1.74 x 1077/2 2.00000 
KM 8 618x107 6.37x 10-7 6.76 x 107-°4 1.74 x 10-7! 2.00000 
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Table 4: Numerical results for fo(«) 
Methods t jer_1| let| Jez+4| |f(v441)| COC 
TM-1 (y=0) 7 3.54x 1071 2.36x 107%! 1.05 x 107 = 1.35 x 107742 2.00000 
(y=- 7 280x107 1.09x 107% 163x107 4.24 x 107-778 2.00000 
TM-2 (y=0) 7 6.03x107!" 6.14x 107%? 6.36 x 1078 = 1.44 x 1077°° 2.00000 
=- 7 248x107 7.30x 10735 6.33 x 107% 6.98 x 107798 2.00000 
Y 
TM-3 (y=0) 7 487x107 3.411079 1.67 x 1078 = 4.95 x 10-?"4 2.00000 
=- 7 640x 10-7! 3.84x 10-4! 1.38 x 10-8! 9.92 x 10-83 2.00000 
Y 
TM-4 (y=0) 7 4.08x 10719 2.05 x 10797 5.15 x 107% 3.29 x 10779 2.00000 
(y=- 7 1.27x 1077?) 1.17x 107-44 1.01 x 10-88 1.68 x 10751 2.00000 
TM-5 (y=0) 7 4.03x 10719 2.01 x 10787 4.99x 1077 2.93 x 107792 2.00000 
(y=- 7 1.25x 1077?) 1.15x 10-44 9.74 x 10789 1.51 x 10751 2.00000 
MNM 7 536x107! 5.56x 10-8! 5.99 x 10-6! 1.49 x 10-739 2.00000 
KM 7 494x107! 4.69x 10-8! 4.23 x 10-® 3.66 x 10-749 2.00000 
Table 5: Numerical results for f3(x) 
Methods t jer-1| let ler+1| |f (t41)| COC 
TM-1 (y=0) 7 849x107“ 6.77x 10-27 4.30 10-53 ~— 8.02 x 10-88 ~—-2.00000 
y=-1 7 205x107 3.46 10-9 9.93 x 10-8 9.12 x 10-18 2.00000 
TM-2 (y=0) 7 2.23x107@ 4.45x10-?8 1.77x 10-5> 6.06 x 10-8! —-2.00000 
y=-1 7 285x107! 6.37 x 107°? 3.17 x 107 = 6.68 x 10-174 2.00000 
TM-3 (y=0) 7 1.91x 107 1.64x 10-36 1.21 x 10-7 1.30 x 107179? 2.00000 
y=-1 7 356x107!" 4.28 x 107°4 6.21 x 107% = 6.26 x 107179 2.00000 
TM-4 (y=0) 7 2.23 x 10-73 3.07 x 10-46 5.81 x 10-82 4.14 x 10-1689 2.00000 
y=-1 7 3.83x 10-79 7.40x 10-49 2.77x 10-7 1.10 x 10-4! 2.00000 
TM-5 (y=0) 7 3.45x 1072! 7.58 x 10-42 3.66 x 10-88 = 1.42 x 10-48 2.00000 
y=-1 7 3.19x 10-2! 5.36x 10-42 1.52 x 10-88 = 3.18 x 10-488 ~—-2.00000 
MNM 7 1.16x10-8 1.27x 10-76 1.53 x 10-52 7.46 x 10-978 ~—-2.00000 
KM 8 7.641077 5.54x 10749 2.91 x 107°? 7.99 x 10-1738 2.00000 


4.1 Unknown multiplicity 


Now, it is clear that the proposed schemes are able to find efficiently the 


multiple roots of nonlinear equations with a known multiplicity. Then, the 


question about unknown multiplicity arises: Is it possible to modify these 


methods in order to simultaneously estimate the roots and their multiplic- 


ity? McNamee [12] presented different techniques due to several authors 


that estimated with unequal results the searched multiplicity. We propose a 
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Table 6: Numerical results for f4(x) 

Methods t er-1| let ler41| |f (t141)| COC 

TM-1 (y=0) 7 3.38x107'7 1.14x107%% 1.30 x 107®° 1.36 x 107®5®5 2.00000 
y=-1 7 252x107" 6.22 10-34 3.78 x 10-87 1.25 x 10-8619 2.00000 
TM-2 (y=0) 7 3.01x107'7 8.96 x 10794 7.95 x 1076” 4.56 x 107®°87 2.00000 
y=-1 7 224x107 4.86 x 10-34 2.29 x 10-7 1.42 x 10-8641 2.00000 
TM-3 (y=0) 7 1.20x 107?! 7.16 x 10-48 2.57 x 10-8 5.21 x 10-845! 2.00000 
y= 1 7 644x 10-72 1.99x 10-48 1.90 x 10-8 6.21 x 10-®5®5 2.00000 
TM-4 (y=0) 7 3.80x10- 1.24x 107% 1.31 x 10-7 2.63 x 10-6988 2.00000 
y=-1 7 269x107 6.06 x 10-86 3.08x 10-7! 7.54 x 10-782 2.00000 
TM-5 (y=0) 7 9.27x 107" 6.70x 10-37 3.51 x 10-7 8.88 x 10-7278 2.00000 
y=-1 7 632x107! 3.04 10797 7.01 x 1077 3.04 x 1077798 2.00000 
MNM 7 348x107 1.21 x 10-33) 1.47 x 10-8 = 2.98 x 10-859 2.00000 
KM 7 348x107!" 1.21 x 10783 1.47 x 1078 =. 2.98 x 1076°6° 2.00000 


modification of our schemes based on the idea of Traub (see [19]), 


m= 


In| f(a)| 


are: f(we) |’ 
— +00 Lt 
Dy | ead 


for estimating the multiplicity in the Schéder scheme 


Tt4+1 = % —™M 


f(t) 
f(z)” 


Expanding this idea to our class of iterative methods, we estimate the mul- 


tiplicity m in expression (5) as 


starting with an initial mp = 


In| f(z2)I 


In [©] 


This provides different expressions 


of m, for each special case of the presented family. These estimators of the 


multiplicity appear described in Table 7. For known methods, MNM and 


KM, the estimation of the multiplicity has been made by the original idea of 


Traub, 
In| f(x 
In Flne,ve] 
and 
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In| f (zz)| 


(1—a) f (m)+af (x) |’ 
In flne.xe] 


Gah Dak 


t= 


respectively. 


Table 7: Estimators of the multiplicity for the special cases of class (5) 
Method mm 
TM-1 0;(1 + aO;) 
TM2 2s 
TM-3 gaethos 
TM4 Gaymnrey? 
ms Seay 


We have made several tests on the same nonlinear equations described in 
Table 2, with the same initial estimations or other more distant estimates, 
using the same stopping criterium and number of digits. The obtained re- 
sults are presented in Tables 8-11. In them, the estimated order of conver- 
gence, COC, does not appear as convergence is linear in all cases for proposed 
and existing methods. We consider that it is due to the natural instability 
provided by divided differences when the estimation of the multiplicity in- 
volves more uncertainty in data. We show the last error estimations, |ez+1|, 
| f(e41)|, the number of iterations, iter, and the last approximation of the 
multiplicity provided by the iterative scheme, m:+1. 

Tables 8-11 show how the number of iterations has increased for all meth- 
ods, possibly due to the instability generated by the estimation of the root 
multiplicity. All the schemes converge with a similar number of iterations, 
and all of them achieve, in all cases, a reasonable estimate of the multiplicity. 
Let us remark that the precision of the estimation of the multiplicity in the 
cases fi(x) to f3(a) is excellent. In the case of f4(a), the relative error of 
this estimation is around 4%. 

Finally, we test our methods with a more challenging function fs(a) = 


4x” +8 sin (x)—47x2+7°—8, whose multiplicity is not obvious. Moreover, there 


7 


exists a wide neighborhood around the multiple roots 5, where the function 


is almost null. Due to the complexity of the problem, the stopping criterium 
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Table 8: Numerical results for f1(2), with zo = 1.8 and estimation of the multiplicity 


Methods iter lez+1| | f (ve41)| Me41 
TM-1(y=0) 32 248x105! 1.49x 10-107 2.0200 
(y = -1) 32. 4.89x 10-5! 5.84 x 107107 = 2.0175 
TM-2(y=0) 32 3.25x 10-5! 2.56x 10-107 2.0174 
(y=-1) 2 s6oxi0-" Waser 2.0175 
TM-3(y=0) 32 4.49x10-5! 4.91 x 10-297 2.0175 
(y= -1) a2  Oriseio= Linx io-** 2.0175 
TM-4(y=0) 32 5.24x 10-5! 6.69 x 10-196 2.0175 
(y=-1) aa 21x10 1e x Or 20175 
TM-5 (y=0) 32 5.12x 10-5! 6.41 x 107197 2.0175 
(y=-1) 32. 7.1se10- «1.256% 10-"° 2.0175 
MNM 32 230K 10-* 1osx10-"” 2.0174 
KM B2- 231% 10- 120810". 2.0174 


Table 9: Numerical results for fo(x), with xo = 1 and estimation of the multiplicity 


Methods iter ler+4| lf (ve41)| Met 
TM-1(y=0) 45 1.161075! 1.79x 10714 1.9605 
(y=-1) 45 1.76x 107°! 4.17x 1071" = 1.9603 
TM-2(y=0) 45 1.42x 10-5! 2.72 x 10-104 1.9604 
(y = -1) 45 215% 10-9 627 x10" 1.9603 
TM-3(y=0) 45 1.351075! 2.47 x 10-104 1.9604 
(y=-1) 45 2.05x 107°! 5.71 x 1071* 1.9603 
TM-4(y=0) 45 1.63x 10-5! 3.59 x 10-14 1.9604 
(y = -1) 45 2.47x 107-5! 8.27x 10-1 1.9602 
TM-5 (y=0) 45 1.80x 107-5! 4.38 x 10-14 1.9603 
(y = -1) 45 2.72x 10-4 1.01 x 10-1 1.9602 
MNM 45 1.15x 10-5! 1.78 x 10-194 1.9605 
KM 45 1.94x 10-5! 5.11 x 10714 1.9603 


have changed to |f(a:41)|< 107°. The numerical results are presented in 
Tables 12 and 13. 
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Table 10: Numerical results for f3(x), with zo = 0.8 and estimation of the multiplicity 


Methods iter lez+1| | f (ve41)| Mte41 
TM-1(y=0) 66 1.01x 10-5! 1.92 x 10-467 8.7317 
(y=—-1) 67 46810- 183x107 8.7333 
TM-2(y=0) 66 3.77x 10-5! 3.07 x 10-482 8.7287 
(y = -1) 67 177K 10-* 324% 10 8.7304 
TM-3 (y=0) 66 1.92x 10-5! 6.68 x 1074% 8.7302 
(y= -1) 67 8.92x 10-5? 6.40 x 10-468 8.7319 
TM-4(y=0) 66 4.75x 10-5! 2.51 x 10-48! 8.7282 
(y=-1) 67 2.22x 107°! 2.58 x 10~*%4 8.7299 
TM-5 (y=0) 67 6.41 10-52 3.19 x 10-469 8.7327 
(y=-1) 67 9.67x 10-5! 1.60 x 10-458 8.7266 
MNM 66 9.79 10->? 1.50 x 10-46" 8.7317 
KM 67 70x10 213%10-7 8.7322 


Methods iter lez+4| | f (ve41)| Mi+1 

TM-1(y=0) 72 3.09x 107°! 2.07 x 107-7579 47.9690 
(y= —1) 72 497x105! 4.75 x 10-7569 47.9610 
TM-2(y=0) 72 3.90x10-5! 2.48 x 10-75% 47.9650 
(y=-1) 72 6.26 x 10-5! 5.57 x 1077555 47.9570 
TM-3(y=0) 72 7.98x10-5! 1.15 x 10-7549 47.9530 
(y¥=-1) 73 5.40 x 10-52. 1.17 x 1077608 47.9980 
TM-4(y=0) 72 7.46x 10-5! 3.98 x 10-755! 47.9540 
(y =-1) 73 5.05 x 10-52 4.06 x 10-7619 47.9999 
TM-5 (y=0) 72 7.16x 10-5! 4.96 x 10-7552 47.9555 
(y=-1) 73 4.84 10-52 5.02 x 10-7611 47.9999 
MNM 72 3.08x 10-5! 1.63 x 107757 47.9699 
KM 72  3.09x 1075! 1.96 x 107757 47.9699 


Table 11: Numerical results for f4(x), with zo = 1.5 and estimation of the multiplicity 


All the methods perform in a similar way, but schemes TM-1 to TM-5 
get better results than MNM and KM methods with xp = 1.8, being approx- 


imately equal in the case where xp = 1. In all cases, we conclude that the 
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Table 12: Numerical results for f5(x), with zo = 1.8 and estimation of the multiplicity 


Methods iter \f(tr41)| Me41 
TM-1 (y=0) 5 5.27 x 10716 3.4652 
(y =-1) 5 6.03 x 10718 3.4857 
TM-2(y=0) 5 5.78 10716 3.4755 
(y=-1) 5 6.13 x 107! 3.4939 
TM-3 (y=0) 5 6.11 x 10-16 3.4910 
(y=-1) 5 6.15 x 10-16 3.4951 
TM-4 (y=0) 5 6.14x 10-16 3.4952 
(y =-1) 5 6.08 x 107'® — 3.4839 
TM-5 (y=0) 5 6.13x 10716 3.4939 
(y=-1) 5 6.12 x 10-'8 3.4909 
MNM 5 5.25 x 10-16 3.4649 
KM 5) 5.26 x 10-16 3.4650 


Table 13: Numerical results for f5(x), with zo = 1 and estimation of the multiplicity 


Methods iter \f(te41)| Me41 
TM-1(y=0) 7 6.05x10-!6 3.4804 
(y=-1) i 5.65 x 10-16 3.4722 
TM-2 (y=0) 7 6.14.x 10718 3.4950 
(y=-1) i 2.18 x 10-'8 3.4440 
TM-3(y=0) 7 5.64x 10-16 3.4719 
(y =-1) 8 5.93 x 10716 3.4635 
TM-4 (y=0) 7 2.86 x 10-16 3.4468 
(y =-1) 8  6.16x 10-6 3.4979 
TM-5 (y=0) 7 1.35x 10716 3.4411 
(y=-1) 8 6.16 x 107! 3.4970 
MNM 7 6.04 x 10-16 3.4788 
KM 7 6.06 x 10-1® 3.4870 


multiplicity of the root is 4. Regarding the precision of the estimation of 


multiplicity, the relative error is, in this case, between 12% and 13%. 
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5 Conclusions 


In this paper, a novel second-order scheme free from derivatives for locating 
the multiple roots of univariate nonlinear functions was developed. A weight 
function L(©;) was used during the development of the new scheme. The 
greatest benefit of using weight functions is that we can produce as many 
special cases as possible and various choices of the weight functions subject 
to some conditions. The error equation in the main result corroborated the 
second-order convergence of the proposed scheme. It has been shown that 
the new family of methods is a generalized one as it contains almost all the 
existing methods of the same type as its special cases. Not only the basins 
of attraction, but also a deep stability analysis in the Riemann sphere have 
also been performed for the new schemes, showing the good performance of 
all the proposed subclasses for small values of the parameter, especially the 
good performance of the TM-2 class. On the other hand, we have assigned to 
each method, parallel to the iterative process to obtain the root, an iterative 
scheme to approximate the multiplicity when it is not known. 

The numerical results for different examples showed that our new method 
is efficient, especially in those cases where the derivative of the nonlinear 
function is not available. 
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